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Abstract 

In this paper we point out the existence of a remarkable nonlocal transformation between the 
damped harmonic oscillator and a modified Emden type nonlinear oscillator equation with linear 
forcing, x + axx + f3x^ + 7a; = 0, which preserves the form of the time independent integral, con- 
servative Hamiltonian and the equation of motion. Generalizing this transformation we prove the 
existence of non-standard conservative Hamiltonian structure for a general class of damped non- 
linear oscillators including Lienard type systems. Further, using the above Hamiltonian structure 
for a specific example namely the generalized modified Emden equation x + ax'^x + (3x'^'^^^ = 0, 
where a, P and q are arbitrary parameters, the general solution is obtained through appropriate 
canonical transformations. We also present the conservative Hamiltonian structure of the damped 
Mathews-Lakshmanan oscillator equation. The associated Lagrangian description for all the above 
systems is also briefly discussed. 
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I. INTRODUCTION 



Dissipative systems are dynamical systems whose phase space volume decreases/ varies 
when the dynamical system evolves in time. This should be contrasted with the conservative 
systems whose phase space volume remains a constant. A dynamical system, which is 
represented by a second order ordinary differential equation (ODE) is said to be a dissipative 
one if the flow function of the equivalent system of first order ODEs turns out to be a negative 
constant^'^. In the case when the flow function becomes zero the underlying system is by 
definition a conservative one^'^. A typical example is the damped harmonic oscillator (DHO) 
equation, x + ax + Xx — 0, where overdot denotes differentiation with respect to time and 
a and A are arbitrary parameters. Working out the fiow function A for this system one 
finds that A = ^ + ^ = —a, where /i and /2 are defined by a; = y = fi{x,y); y = 
—ay — Xx = f2{x,y), which confirms that the system under consideration is a dissipative 
one. On the other hand, for a conservative system, we have a set of first order equations, 
X — ^ — fi{x,p)] p — — f2{x,p), where H is the Hamiltonian function, and the fiow 
function is always equal to zero which can be straightforwardly confirmed by substituting the 
Hamilton equations of motion. This observation reveals the fact that if one is able to find a 
time independent integral for a dissipative system then this time independent integral can be 
correlated to a Hamiltonian function, perhaps in a non-standard form in terms of a new set 
of canonical variables, for the second order equation under consideration. Such a conserved 
Hamiltonian description also leads one to further investigations, including quantization^. 

Since the pioneering work of Bateman^ the quest for a Hamiltonian description for the 
damped harmonic oscillator equation was pursued by several authors and only very recently 
conserved Hamiltonian description for this system^ for all the three parametric regimes have 
been constructed by three of us . A possible quantization has also been suggested. 

In a parallel investigation we have also found that another dissipative type system, namely 
the modified Emden type equation (MEE), also admits time independent integral and a 
conservative Hamiltonian description^. A question which now naturally arises is whether 
the aforementioned systems are isolated examples or there exists a wider class of dissipative 
systems that admit conservative Hamiltonian description. If so, how to isolate and classify 
them? A more detailed investigation in this direction reveals the fascinating fact that one 
can map the damped harmonic oscillator equation onto the modified Emden type equation 
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through a nonlocal transformation. Interestingly the time independent integrals and the 
Hamiltonian description for the nonlinear system can also be derived from the linear one by 
simply substituting the same nonlocal transformation at the appropriate places. 

In this paper, after presenting the above novel results, we report a rather general transfor- 
mation which can map the damped harmonic oscillator to a larger class of damped nonlinear 
systems that admit conservative Hamiltonian description. We present the general nonlinear 
oscillator equation and explicit forms of the Hamiltonian. The Hamiltonian forms obtained 
by this procedure are of non-standard (not equal to the standard potential plus kinetic en- 
ergy) forms. Recently considerable interest has been shown in the classification nonlinear 
dissipative equations which admit non-standard conservative Hamiltonian structure 
We then consider a specific nonlinear system, x + ax'^x + /Sx^^"*"^ = 0, where a and (5 
are arbitrary parameters, which is a natural generalization of the DHO equation and the 
MEE and discuss the dynamics/integrabihty of this equation in some detail. Explicit so- 
lution for this equation was constructed for the specific parametric choices^^'^^ (5 = j^^^ 
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and /3 = 4(^5:1) • For the former choice of parameter we have recently shown that this 
equation can be linearized to the free particle equation through a generalized linearizing 
transformation^^. Also the q = 1 case (MEE) has been completely integrated by using the 
underlying Hamiltonian structure^ for all values of the parameter a and (3. However, neither 
the Hamiltonian structure nor explicit solutions for arbitrary choices of g, a and j3 for this 
equation have been reported and we present the results here. Finally we also present the 
Hamiltonian structure of the damped Mathews- Lakshmanan oscillator, whose undamped 
version exhibits amplitude-dependent harmonic type oscillatory solutions''^^'''^^. 

The plan of the paper is as follows. In section II we discuss the time independent integrals 
and Hamiltonians of the two dissipative systems, namely the DHO and MEE with linear 
forcing and deduce the nonlocal transformation which interrelates each other. In section HI 
we introduce a more general nonlocal transformation and substitute it into the DHO equation 
and construct the nonlinear generalizations of the DHO. We also derive the associated time 
independent integrals and Hamiltonian for this general nonlinear oscillator equation. In 
section IV, we consider as a specific example, namely the generalized MEE, and discuss the 
Hamiltonian structure and obtain its general solution by integrating the canonical equations 
of motion after suitable canonical transformations. In section V, we present as a second 
example the time independent integral of motion and conservative Hamiltonian structure of 
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the damped Mathews-Lakshmanan oscillator. In section VI, we briefly discuss the associated 
Lagrangian description for all the above systems. Finally, we present our conclusion in 
section VII. 



II. HAMILTONIAN STRUCTURE OF DHO AND MEE 

In this section we briefly recall the Hamiltonian dynamics associated with the two dis- 
sipative systems, namely the damped harmonic oscillator and the modifled Emden type 
equation, and then show how they are interrelated. Note that the former one is a linear sys- 
tem and latter one is a nonlinear system. One may essentially consider the transformation 
to be introduced as a linearizing transformation of the latter. 



A. Damped harmonic oscillator 



To start with let us consider the damped harmonic oscillator (using a different notation 
for convenience of comparison) 



y" + ay' + Ay = 0, 



d_ 



(1) 



where a and A are arbitrary parameters. Recently, we have identifled the following time 
independent integral of motion^ for the system ([T]) : 

' (r - 1) / , a \ , , (r - 1) , 

7 — ^ iy' + -y) {y' + -(^y) 

[r — 2) \ r J r 



(l-r) 



^ log[2/'2 + ayy' + Xy^] + ^ tan"^ 



a^ > 4A 



ay' + 2Xy 
2uy' 



a2 = 4A 



, a^ < 4A, 



(2) 



where r = ^(ct ± ^ c? — 4A) and uj = |V4A — a^, for the overdamped, critically damped 
and the underdamped oscillations, respectively. From these time independent integrals we 
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have also derived the following time independent Hamiltonian ^, 



2) 



(r-2) 
(r-1) _ 



(r-1] 



apy, p > 0, > 4A 



^ - { log(p) - ^apy, 



p> 0, a'^ 



4A 



(3) 



^\og[y^sec^{upy)] - ^py, 
where the canonical conjugate momentum p is defined by 



p 



y + 



ir-l] 



-ay 



(l-r) 



1 



tan 



V + ay 
2ujy 



a^ < 4A 



a^ > 4A 



a^ < 4A 



(4) 



From the Hamiltonian ©-dl]) one can straightforwardly write down the canonical equations 
of motion for all the three regimes and obtain the known exact solutions straightforwardly. 
Note that no multivaluedness arises in the system due to the constraints on the momentum. 
Case 1,2 > 4^ 



Substituting the first of the expressions in ([3]) into the canonical equations y' 



dH ^1 



dp 



p = - 



m 

dy ' 



one gets the following equivalent system of first order ordinary differential equations for the 
damped harmonic oscillator in the overdamped and critically damped {a^ > 4A and = 4A) 
parametric regimes, 



y' = (p) 



(r-1) 



-ay, p 



1) 



-ap. 



(5) 



One can easily check that the second order equivalence of this equation coincides exactly 
with ([T]) and that the standard form of the solutions for the overdamped and the critically 
damped cases follow naturally by integrating the above system of first order ODEs. We 
mention here that the divergence of the fiow function of the underlying canonical equations 
of motion is zero, A = ^ + ^ = 0, which in turn confirms that the damped harmonic 
oscillator has a non-standard conservative Hamiltonian description. 
Case 3 < 4A 

Finally, for the underdamped case, one can write down the Hamilton equations of motion 
in the form 



y' = ^l/(2u;tan(wp?y) - a), p = ^p{a - 2utan{upy)) - ^. 



(6) 



Again rewriting this equation as a single second order equation in the variable y one ends 
up in the form ([1]). The divergence of the flow function in this case is also zero. 

The existence of time independent Hamiltonians ([3]) in a non-standard form and the 
basic definition of the divergence of flow function lead us to the conclusion that the damped 
harmonic oscillator has a non-standard conservative Hamiltonian description as well, in spite 
of its known dissipative nature. 



B. Modified Emden type equation 

In a parallel investigation we have also found that the MEE, x + axx + f3x^ = 0, admits 
a time independent conservative Hamiltonian description^. The MEE, which is also known 
as the Painleve - Ince equation, is an extensively studied and physically significant equation 
in the contemporary nonlinear dynamics literature^^"^^. On the other hand the MEE with 
linear external forcing, 

X + axx + Px^ + 'jx = 0, (7) 
possesses certain unusual nonlinear dynamical properties for the parametric choice (3 = 
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We have shown that this system admits nonisolated periodic orbits of conservative 
Hamiltonian type^. These periodic orbits exhibit the unexpected property that the frequency 
of oscillations is completely independent of amplitude and continues to remain as that of 
the linear harmonic oscillator^. 

The more general system ([71) with a, (3 and 7 being arbitrary parameters also admits the 
following time independent integrals of motion^^, 

x + ^x^+ j 1 ( X + ^^^x^ + '"^\ a- > 8/3 



r — 2) \ 2r a{r — 1) J \ 2r a 



^ \og[x + ]a{x'^ + ^)], a^ = Sp 



I = < 



{x + \a{x^ + ^)) 4 ^ a^. 



hog[i^ + ^x{x^ + l) + lix^ + lr] 



(8) 



ax + 2f3{x^ + ^) 



2ujx 



a^ < 8/3, 
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where r 



4/3 



a 



± — 8/9) and io = | a/8/3 — a^. From the time independent integrals we 



have identified the following Hamiltonian for the system ([7j), 



H 



ir-1) . 



(r-2) 



2r 



2 P^l 

apx , 

a 



^ ^ 47 



X 



a 



a" 



a 



(3' ^ 2 

where the canonical conjugate momentum p is defined by 



2^ 
P' 



p 



X + 



1^ 



2r 



2 , n 

ax H 

a 



(l-r) 



ujix"^ + i 



tan ""^ 



4x + a(x2 + 
2a;(x2 + ^) 



The Hamilton equations of motion follows straightaway from (jH]) as 



Cases 1 & 2 x 
Case 3 x ■ 

p 



pl-r 
1 



2r 



— , P 



a(r — l)px 



4/3 



(/?x^ + 7) 2u; tan 



1 / 2 , 7^ 



r 

— a 



x{pa{j3x'^ + 7) - 4/3 - 2puj{j3x'^ + 7) tan[|pu;(x^ + ^)]) 



(9) 



(10) 



(11) 



(12) 



2(/3x2 + 7) 

Note that in all the above three cases, the divergence of the fiow function A is zero, corre- 
sponding to a a conservative Hamiltonian description. However, unlike the damped harmonic 
oscillator one cannot integrate the canonical equations (fTTIl and ( fT2l) and obtain the solu- 
tions straightforwardly. To overcome this difficulty one should introduce suitable canonical 
transformations and change the Hamilton equations flTTl) and f[T^ into simpler forms so that 
they can be integrated in a straightforward manner as demonstrated in the case of MEE 
(with 7 = 0) recently by us^. 



C. Transformation connecting DHO and MEE 

By comparing the structure of the integral of motion, Hamiltonian function and canonical 
equations of the damped harmonic oscillator with that of the equation ([7]), namely equation 
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with ©-dH) with ©-([lO]) and ([S])-© with (in])-(IT2]), one can identify that the two 
systems are transformed into each other through the nonlocal transformation 

y = y + -, dr = xdt, (13) 

with the identification A = 2/3. One can also check directly that the above transformation 
also maps the damped harmonic oscillator equation ([T]) onto the MEE with forcing, equa- 
tion (^^, and vice versa. Consequently one can treat the transformation as a linearizing 
transformation (albeit nonlocal) of the nonlinear equation ([7]). 

Note that the transformation (fT3l) is not the only possible linearizing transformation at 
least for specific parametric choices. We mention here that one can also transform the 
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nonlinear system ([71) with P = ^ into a linear harmonic oscillator equation {U + jU = 0) 
by introducing another nonlocal transformation of the form ^'^^ 

jj ^ ^J^axii)di^ r = t. (14) 

The nonlocal transformation ffl^ is different from (fT3!) in the following respect. In (fT3|) 
the nonlocality is introduced in the independent variable whereas in f|T^ the nonlocality is 
introduced in the dependent variable. We may also add that the nonlocal transformation 
f|T3|) has some similarity with the well known Kustaanheimo-Stiefel transformation used 
in atomic physics^'^. Even though both the nonlocal transformations fll3p and (1141) map 
the nonlinear equation into a linear one and vice versa, the nonlocal transformation of the 
type f|T3|) is much useful in identifying Hamiltonian structures associated with the nonlinear 
system whereas the nonlocal transformation of the type (fT^ is more useful in constructing 
general solution for the transformed nonlinear system. For more details about the nonlocal 
transformation of the type (1141) we may refer to Ref. 22. In the following we confine our 
attention to the nonlocal transformation of the form f[T51) and its generalization. 



III. A GENERAL CLASS OF NONLINEAR DAMPED OSCILLATOR : HAMIL- 
TONIAN DESCRIPTION 

We find that the transformation f|T3|) is a specific case of a rather general nonlocal trans- 
formation of the form 

y = [ f(x)dx, dr = ^\ \ dt. (15) 
J 9{x) 



For example, restricting f{x) = x and g{x) = 1 in (fT5|) one gets exactly (fT3|) . 

The nonlocal transformation f|T5l) modifies the damped harmonic oscillator equation ([1 
to the general class of nonlinear oscillators of the form, 

9'{x) .2 , f{x) . f{x) 



X + 



-X + a——-x + A- 



f{x)dx = 0, 



(16) 



g[x) ' g{x) g{xY 

where /(x), g{x) are arbitrary function of x. 

Applying the above nonlocal transformation ( fT5l) to the damped harmonic oscillator equa- 
tion ([T]) we obtain the following time-independent integral of motion for the nonlinear system 
(fT6|). that is. 
Case 1. > 4A 



/ = ^ lg{x)x H — J f{x)dx\ lg{x)x H a / f{x)dx 



(l-r) 



Case 2. ^2 = 4^ 



+ |a / f{x)dx 



^og[g{x)x + y /(a;)'^^;] 



Case 3. < 4A 



log[(7(a;) i + ag{x)x / f{x)dx + A( / f{x)dx) ] 



a _, 

H tan 

2cj 



a(7(x)a; + 2A / f{x)dx 
2ujg{x)x 



(17) 



(18) 



(19) 



where u> = |V4A — and r = ^(a ± \J — 4A). 

The Hamiltonian for the equation ( |T6l) can also be constructed by simply substituting 
the transformation ffTSl) into (131) so that the latter reads 



f r — 2 

(r — 1) / P \ '■"^ a('" — 1) P 



(r-2) \g{x) 



r g[x) 



f{x)dx, 



H = < 



a 



2g{x 



-p J f{x)dx + log 



I 2 



log 



f{x)dxY sec^ 



9[x) 
P 



up 
g{x) 



> 4A 



= 4A 



(20) 



f{x)dx 



ap 
2^) 



f{x)dx, < 4A, 
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where p is the canonical conjugate momentum defined by 



g{x)x+ J f{x)d. 



n r-1 ' 



X 



> 4A 



(21) 



tan 



2g{x)x + a f f{x)dx 
2uj J f{x)dx 



< AX. 



One can easily check that the canonical equations of motion take the form 



a' > 4A 



X 



p 



9{x) 



a / f{x)dx 



P 



{r - l) ^^ f{x) _^ pg'{x 



< 4A 



X 



r " g{x 
J f{x)dx 



2g{x) 



2uj tan 



g{xf 
pu 



P 



a{r — 1) 



g{x 

f{x)dx 



P 



2g{xY J f{x)dx 
xp / f{x)dx J + p 



9[x) 
f{x)g{x)[ 2g{x) 



a 



a — 2uj tan 



f{x)dx 



pu J f{x)dx 



(22) 



f{x)dx 



a — 2a; tan 



9{x] 
pu J f{x)dx 
g{x) 



g'{x) 



(23) 



One can straightforwardly check that the second order equivalence of the equations (!22|) and 
f l23|) coincides exactly with f|T6|) in the appropriate parametric regimes. Further, we find that 
the fiow function for the canonical equations of motion is zero in all the three parametric 
regions, equations ( I22|) -(l23l). which in turn confirms that ( |T6i) or ( l22l) -( l23|) has a Hamiltonian 
description. Further, the integrability of f|T6|) is automatically ensured by the existence of 
the time independent Hamiltonian (!20|) in the Liouville sense. 

We also note here that equation f|T6l) can be transformed into an Abel equation of the 
second kind^^, 



/ , 9'{x) 2 , f{x) , . fix) 
WW H r-rW + a— — + A , 



g{x) 



f{x)dx = 0, 



(24) 



'g{x) g[x) 

through the transformation w{x) = x. Then the time independent first integrals of ([24 
with g{x) = constant can be constructed by the procedure described in Ref. 25. In the 
case where J f{x)dx is an invertible function of x, then the procedure given in Ref. 25 
to construct the solution will lead to quadratures in terms of certain complicated integrals 
which however cannot be evaluated in general. For further details one may refer to Ref. 26. 
Therefore one should adopt a different procedure to construct the general solution of (fT6|) . 
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The effective way to proceed further is to transform the Hamiltonian (120|) to a simpler form 
through suitable canonical transformations as demonstrated recently by us for the case of 
the MEE^. However, to adopt this procedure one should specify the explicit forms of f{x) 
and g{x). We demonstrate the procedure with a specific example in the following section. 



IV. GENERALIZED MEE : HAMILTONIAN STRUCTURE AND GENERAL SO- 
LUTION 

To illustrate the ideas given in section 3, we focus our attention on the case f{x) = x"^, 
g{x) = 1 so that equation (fT6|) now becomes 

X + ax'^x + Px'^'^^^ = 0, (25) 

where (3 = The reason for choosing this form of f{x) for illustration is that it provides 

a natural generalization of the damped harmonic oscillator and the modified Emden type 
equation. 

Substituting f{x) = x'^ in equations ( |T71) - ( fT9|) we obtain the following time independent 
integral of motion for the equation (l25l) . that is 
Case 1. > 4/5(g + i) 

Case 2. = 4/3(g + 1) 



x+fx^+i ' 2 
Case 3. < 4^(g + i) 

/ = - \og\x^ + axx'^^ + + — tan-i 
2 ^ (q + l) 2uj 



- \og[x + ^x^+i] (27) 



ax + 2/32;'?+^ 
2ux 



(28) 



where r = 2/3{g+i) i'^ ^ V'a^ - 4/3(g + 1)), u = \^Aj3{q + I) - and a = The Hamil- 
tonian for the equation (l25l) can be deduced from equation (1201) which turns out to be 
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> 4/5(g + 1) 



H 



log(p) 



a 



a 
—px 



9+1 



-px 



- log 

2 ^ 



a; 



= 4p{q + 1) 



(29) 



where the canonically conjugate momentum is defined by 



P 



1 



(r-1) 
('Z + l) 



a; + 



(r-1) 



-ax 



9+1 



(1-r) 



(30) 



g+i 



tan ^ 



ax 



''+1 + 2(g + l)x' 



2a;x5+i 



Here we note that the integrals fl26l) - fl28l) can also be derived systematically through 
various methods available in the recent literature. To name a few, we cite Prelle-Singer 
procedure^^, method of reducing Lienard equation to the Abel equation form^^, factoriza- 
tion method^^'^^ and so on. However, all these methods provide only time independent 
integrals but neither the integrals nor the Hamilton canonical equations can be integrated 
straightforwardly. 



A. Method of obtaining general solution 

In this section we transform the Hamiltonian fl29|) into a simpler form, by introducing a 
suitable canonical transformation and then obtain the solutions. 

Case 1: a^ > 4/?(g + 1) 

The Hamiltonian for this parametric regime is given by 

(r-1) (L^ (r-1) ^ „ , 1 (a^ ± aJa'^-4:P(q + l)) „ 

By introducing a canonical transformation 

x = g(^^y, p=l{U'^~'P^^'y (31) 
12 



the Hamiltonian (l29l) can be transformed into the form 
where we have defined ai = {^r-^ (|) '"^"^^ ) '^i = ■^^^^^rrij^) 



(32) 



= fa and 



The canonical equations of motion for the transformed Hamiltonian ( ]32l) now reads 

f/ = ainit/^^F^i-S (33a) 



"1 rrmi-l 



27/1 f/). 



(33b) 



One may observe that for the choice ni = 1 equation fl33|) becomes uncoupled and one can 
integrate the resultant equations straightforwardly. In the following, first we consider the 
general case and then discuss the particular case ni = 1 which corresponds to the parametric 
choice P 



(9+2)2 



Substituting fl33al) into (!32|) and rewriting the latter, we get 



E = ^ 
(cTirii)"!-! 

Rewriting equation fl34|) we obtain 



+ ViU. 



(34) 



U 



■[E-rj,U^ 



ni ^ 1 



(35) 



where ai 



ni ■ Integrating the above equation one obtains. 



t-tn 



(36) 



where m = and = The above integral can be split into two cases depending 

upon the values of m and ra, that is 

Ijm-l 



r]i{2n - m - 1)(^ - r^i?/^)""! 

(m- 1)^ 



?7i(2n — m — 1) J {E — rjiU"^) 



m-2 



-dU 



m, n > 



f - to = < 



(37) 



1 



E{1 - m)U>^-^{E - ?7i[/2)"-i 

^ ?7i(m + 2n — 3) 



£;(i-m) y U'^-^iE -riiU^y 



-dU m < 0. 
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The integrals on the right hand sides of (1371) can be integrated again and again until all the 
integrals are exhausted, thereby giving the general solution in an implicit form. 

Now we analyze the case ni = 1 in (155]) . For this choice, equation (I33al) can be integrated 
to yield U in the form 

?7 = [(l-mi)(ait + ci)]^, (38) 

where Ci is an integration constant. From equation ( l32i) one can express P in terms of U 
and E and substituting the latter into the first expression in fl3T]) one gets 

(39) 



where ai = ^ (I) ' ' ~ ~ 0^ 2(9+2) ^"^^ "^1 ~ f+T' Substituting fl38|) into fl39|) and 
simplifying the resultant expression we arrive at 



x(t) = I 2— g^"^(g + iHg + 2)M + cO^ , 

^ ^ ' 2^(g + 2) + a2«+ig«+i(ait + Ci)''+i ' ' ^ ^ 

with E and ci are two arbitrary constants. We mention here that the linearizable case 
= 9/3 belongs to this case and the solution exactly agrees with the known one in the 
literature^°~^^. 

Case 2 : a^ = 4/3(g + 1) 

Using the same canonical transformation, equation fl3Tl) . one can transform the Hamilto- 
nian, H = log(p) — —px'^'^^, into the form 

H = log(?7) + ^ log[P] + rj2U^ = E (41) 

q q 

where 772 = —^q^''~^^\ The corresponding canonical equations become 

E>-^. (42a) 

P = (^) ^ + ir^a (42b) 

From fl42al) one can express P in terms of U and substituting this into fl4ip one can bring 
the latter to the form 

E = \og{U) - log(f7) + r]2U^ (43) 
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which in turn gives us 



U = exp[ 



qr]2U'^ 
(9 + 1) 



(44) 



where Ei = exp[^^^]. Integrating the above equation, one obtains the solution in the form 



t-to = — 
El 



dU 



2Ei 



1 qr]2lP_ 
g + r q + l 



(45) 



where F is the gamma function^^. For the choice q = 1 the integral (H5l) can be evaluated 
in terms of error function 



Case 3 : < A(3{q + 1) 

Now we focus our attention on the underdamped case. Using the canonical transformation 
( !3T|) . we rewrite the underlying Hamiltonian in the form 



The associated canonical equations read 

q+l 



(g+i) 

' sec[ . " .. [/^ 



^2(g + i; 



u 



p 



qP ' 

2quU' tan[^] - ((g + l)q^+'U^a - 2{q + 1) 
2q{q+l)U 



From fl47ap one can write P = ^ti- and substituting this in the Hamiltonian 
simplifying we get 



(46) 



(47a) 
(47b) 
and 



H = log[[/[/^ sec[ ^ 



2(g+l) 



U^]] + V2U'' = E. 



Rewriting equation (1481) we have 



U 



iexp[^(E-r^2f/' 



U 



sec 



2{g+l) 



f/2] 



9 + 1 



(48) 



(49) 
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Integrating equation (H9i) . we get 



t-tn 



U seel 



2(g+l) 



rdU 



exp[^(i?-r^2t/2)] 

g 0^(0; — 2iqi) (3g + 2)a; — 2iqqi 



g + 1' 2(g+l) 



2(g + l)a; 



-exp 



(cosh[g2] - sinh[g2]), 
where gi = ?72(g + 1), g2 = ^ and F is the hypergeometric function^°. 





sec 


_2{g+l)_ 


(l + exp[i^])^ 


"J 


q{2{q + 1)7]2 + icu) 



(50) 



V. DAMPED MATHEWS-LAKSHMANAN OSCILLATOR 
STRUCTURE 



HAMILTONIAN 



As a second example we consider the following damped version of the Mathews- 
Lakshmanan oscillator equation^^, 



X 



,x^ + 



tX + 



0, 



(51) 



1 + Aix2 1 + Aix2 1 + Aia;2 

where we have included a nonlinear damping term to the Mathews-Lakshmanan oscillator 
The Mathews-Lakshmanan oscillator (equation (ISTl) with a = 0) possess simple trigonomet- 
ric solution^^'-^^, x{t) = Acos^lt, fi^ = ^_^^^^2 with the Hamiltonian H = p'^(l+Xix'^) + Y^^^ 
. We now find that the damped case (IHTj) . with a 7^ is also integrable and admits a con- 
servative Hamiltonian structure. 

Comparing the above equation with f|T6|) . we find 

1 . 1 



The integral of motion is then 



VT+Xlx^ 



fix) 



[1 + Aix2 



(52) 



(r — 1) (rx + ax) 



(r - 2) rVl + Xix' 



rx + ax[r 



2x 



2x — X 



a 



log 



2i; + aa; 



l-r 



a2 >4A 



2V1 + Aix2_ 
+ 2Ax 



4A 



(53) 



2ljx 



+ axx + Ax^ 
(1 + Aix2) 



< AX. 
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Substituting the above forms of g{x) and f{x) in the expression for the Hamiltonian (l20l) . 

we get 



' (r-l) r /^— ; — a{r-V 



H 



(r-2) 



■px, o? > 4A 



a 

< —px + log 



p 



^log 



sec^ [cjpa;] 



where the canonically conjugate momentum 



a^ = 4X 



— x, < 4A, 



(54) 



1 



rx + ax(r — 1] 



V1 + Aia;2 V rVr+Alx2 



l-r 



p = < 



■ tan ^ 



a + 2i; 
2a;x 



1 



a2 > 4A 



< 4A. 



(55) 



Using the canonical transformation P = - log(x^), U = px, one can rewrite the Hamiltonian 
for the underdamped case (a^ < 4A) as 



H = ^ {log[sec\u;U)] - log[e-^^ + Ai] - aU) = E. 



(56) 



Following the same procedure illustrated in the previous example, one can arrive at the 
following equation 



U =1- cos^{uU) exp[2£^ + aU]. 

Obviously this can be rewritten as the quadrature 

dU 

t-tn 



(57) 



(58) 



1 - cos^iioU) exp[2E + aU] ' 
The above integration can be easily performed for the choice a = and one can recover 
the known solution of the Mathews-Lakshmanan oscillator. For a 7^ 0, the integration is 
nontrivial and requires further investigation. Similarly, for the overdamped and the criti- 
cally damped cases, one can use the canonical transformation U = smh~^ (^/X^x) , P = 
pyjl + AiX^ to reduce their corresponding Hamiltonian to simpler forms. However, again we 
are left with a quadrature for the choice a 7^ 0. Work is in progress to find other suitable 
canonical transformations to find the general solution explicitly for the case a 7^ 0. 
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VI. LAGRANGIAN DESCRIPTION 



We now present the equivalent Lagrangian description for the systems studied so far. One 
can follow two different procedures. Since both the Hamiltonian, H , and the corresponding 
canonical momentum, are available for each one of the systems studied in the previous 
sections, one can write down the associated Lagrangian as 



L = px — H. 



(59) 



Alternatively, from the known form of the Lagrangian for the damped harmonic oscillator 
given by us in Ref. 5, and applying the nonlocal transformation f|T5l) . one can deduce the 
corresponding Lagrangian for the nonlinear damped oscillator equations. One can easily 
check that both the methods give rise to the same Lagrangian. 

Considering now the damped harmonic oscillator ([T]) and its Hamiltonian form ([3]) along 
with the conjugate momentum (j4j), the Lagrangian can be expressed as 



L 



y 



log 



2uj 



-ay 



(2-r) 



> 4A 



y' + ^oiy 



a tan ^ 



a 



4A 



(60) 



2uy' 



— tan 

y 



ay + 2j/' 
2ujy 



+ -log [y'^ + ay'y + \y^] 



< 4A. 



Now applying the nonlocal transformation (ITSl) to the above Lagrangian (!60!) . one can 
readily obtain the Lagrangian associated with the general nonlinear damped oscillator equa- 
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tion (dSD as 



g{x)x + 



-a 



f{x)dx 



(2-r) 



> 4A 



log 



g{x)x + -a 



a 

— tan 
2uj 



f{x)dx 
ag{x)x + 2A / f{x)dx 



a 



4A 



(61) 



2ujg{x)x 



g{x)x 



uj J f{x)dx 
1 



tan 



-1 



+ -lo£ 
2 ^ 



a J f\x)dx + 2g{x)x 



2uj J f{x)dx 
g{xYx^ + ag{x)x I f{x)dx + A( / f{x)dxY 



< 4A. 



Naturally, the same form follows from the Hamiltonian fl20|) and the conjugate momentum 
dH]) using (jSni). The Lagrangian for the specific examples considered in Sees. IV and V can 
be deduced by specifying the forms of f{x) and g{x) in equation flM]) . It may be noted that 
all the above Lagrangians are of non-standard type, that is of the forms which cannot be 
written in the standard form as 'kinetic energy' minus 'potential energy'. In particular, one 
can readily check that the non-standard Lagrangian 

1 



i; + I / b{x)dx ' 



(62) 



deduced by Musielak^ for the nonlinear ODE 



X + b{x)x + -b{x) I b{x)dx 
9 



0, 



(63) 



follows from the above general form (I6T!) by choosing f{x) = b{x), g{x) = 1, a = 1 and 
A = I in equation flTBl) . Note that the above equation fp5]) itself includes the MEE discussed 
in Refs. [8,9] as a special case with b{x) = kx. 



VII. CONCLUSION 



In this paper, we have investigated a class of nonlinear dissipative systems which ad- 
mit a non-standard time independent Hamiltonian description through a novel nonlocal 
transformation. The procedure is simple and straightforward. By introducing a nonlocal 
transformation in the 'source equation', namely the DHO equation, we are able to generate 
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a class of 'target equations', namely the nonlinear generalizations of DHO. We have used the 
same nonlocal transformation to deduce the time independent Hamiltonian for the nonlinear 
equation. The nonlocal transformation introduced in this paper is different from the one 
which we have adopted in Ref. 22 and has certain salient features in identifying the time 
independent Hamiltonian structure. To illustrate the procedure we have considered two 
specific systems, namely the generalized MEE and the damped Mathews-Lakshmanan oscil- 
lator equation. We obtained the solution of generahzed MEE by integrating the canonical 
equations of motion after introducing a suitable transformation. The associated Lagrangian 
description for all the above systems is also briefly discussed. A similar analysis can also 
be performed in the case of two degrees of freedom systems. The details will be presented 
separately. 
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